Conditional preparation of states containing a definite number of photons 



Malcolm N. O 'Sullivan, ^'[^jKam Wai Clifford Chan,^ Vasudevan Lafcshminarayanan,^ and Robert W. Boyd^ 

^ The Institute of Optics, University of Rochester, Rochester, NY 14627, USA 
^School of Optometry, Department of Electrical and Computer Engineering, 
and Department of Physics and Astronomy, University of Waterloo, Waterloo, ON, N2L 3G1, Canada 

(Dated: February 4, 2008) 

A technique for conditionally creating single- or multimode photon-number states is analyzed using 
Bayesian theory. We consider the heralded A'^-photon states created from the photons produced by an 
unseeded optical parametric amplifier when the heralding detector is the time-multiplexed photon- 
number-resolving detector recently demonstrated by Fitch, et al. [Phys. Rev. A 68, 0438f4 (2003).] 
and simultaneously by Achilles, et al. [Opt. Lett. 28, 2387 (2003).]. We find that even with 
significant loss in the heralding detector, fields with sub-Poissonian photon-number distributions 
can be created. We also show that heralded multimode fields created using this technique are more 
robust against detector loss than are single-mode fields. 



A Fock state of the electromagnetic field contains a def- 
inite number of photons and hence displays no intensity 
fluctuations, fn contrast, the distribution of the inten- 
sities of a classical field of light has a variance at least 
as large as its mean [T]. Thus, Fock states and other 
states with sub-Poissonian photon-number statistics are 
of interest not only because of their fundamentally non- 
classical nature but also because of their potential to 
make radiometric measurements with greater accuracy 
than possible using classical states of light. Research in 
Fock state generation is also motivated by its applications 
in quantum information science [2] including the areas of 
quantum cryptography [3] and quantum computing 0] . 

Unfortunately, it can be experimentally challenging to 
produce light fields that contain definite numbers of pho- 
tons. Several theoretical proposals aimed at creating 
Fock states "on demand" have been offered [5l|6j|7], but 
they involve trapped atoms in high-Q cavities and to the 
best of our knowledge only single-photon states have been 
created in the laboratory [8J. Several stochastic sources 
of one- and two-photon states have been experimentally 
demonstrated using various processes including atomic 
and molecular fluorescence [TO] , Coulomb blockade for 
electrons [TT] and action in a micromaser [T^]. In addi- 
tion, a particularly simple method of generating single- 
photon states using the photon pairs created in the pro- 
cess of spontaneous parametric down-conversion was first 
demonstrated by Hong and Mandel |f 3j . In their experi- 
ment, they used the single- photon detection of one pho- 
ton of the pair to "herald" the presence of the other pho- 
ton of the pair. This heralding scheme is generalizable 
to create multiphoton states of light. In an unseeded 
optical parametric amplifier (OPA), pump photons in- 
teract nonlinearly with a noncentrosymmetric crystal to 
produce multiple pairs of less energetic photons termed 
the signal and idler through the process of parametric 
down-conversion (PDC). In general, the number of pho- 
tons found in the signal and idler fields is random and 
known to obey thermal statistics, which can be charac- 
terized by the gain g of the OPA. But since the emission 



of every signal photon requires the simultaneous emission 
of an idler photon, if one knows the number of photons 
in the idler field then the number in the signal field is 
determined with complete certainty. Thus, in principle, 
by allowing transmission of the signal field only when 
the idler field is found to contain n photons, the signal 
field will be guaranteed to also contain exactly n pho- 
tons. This idea was considered by Holmes et al. [14] for 
ideal photon-number-resolving detection. 

However, the photon-number-resolving detection 
needed to create heralded multiphoton states is not 
trivial to implement and is an area of active research. 
Photon-number-resolving detection has been achieved us- 
ing superconducting transition-edge sensors |15j . visi- 
ble light counters [TB], charge integration photon de- 
tectors [T7], and by time-multiplexed detection involv- 
ing only commercially available APDs [IHl US]. In the 
present work, we show that, in particular, the photon- 
number resolution achieved using time-multiplexed de- 
tectors (TMDs), as first conceived by Fitch et al. [T5] 
and independently by Achilles et al. [19 , is well-suited 
for creating heralded fields with near-definite numbers of 
photons even in the presence of significant detector loss. 
This technique has been demonstrated experimentally for 
two-photon states [5^ but has not yet been analyzed in a 
comprehensive fashion for states containing greater num- 
bers of photons. 

A schematic of the proposed heralding technique is 
shown in Fig. [l] As mentioned above, the source of 
photon pairs is the output of an unseeded OPA and 
the heralding detector is a TMD. The TMD provides 
photon-number resolution by first randomly splitting the 
incoming state into many pathways that are separated in 
time from one another, i.e. they are separated into time 
bins. The number of photons in the incoming field can 
be estimated to be equal to the total number of photons 
detected by two avalanche photodiodes (APDs). Any 
errors in this estimation will introduce uncertainty into 
the number of photons contained in the heralded state. 
Undercounting is the main source of errors in TMDs. 
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FIG. 1: (a) Scheme to produce heralded fields containing 
a definite number of photons. An unseeded optical para- 
metric amplifier (OPA) produces photon pairs through the 
process of parametric down-conversion (PDC). When a time- 
multiplexed photon- number-resolving detector (shown in (b)) 
detects Ui photons in the idler field, it heralds the presence of 
Us photons in the signal field with probability Psig{ns\ni) . 



The two primary mechanisms for undercounting errors 
are losses and inadequate splitting of the incoming beam. 
Losses in the detection system stemming both from the 
beam splitters and from inefficiencies in detectors them- 
selves can be quite significant and must be included in 
the analysis. Inadequate splitting occurs because of the 
finite number of times the incoming beam is split. Hence, 
a chance exists that two or more photons will occupy the 
same time bin resulting in an undercounting error at the 
APDs, which can perform only on/off detection. More 
quantitatively, a system containing m beam splitters will 
have M = 2™ different time bins. If the number of in- 
coming photons N is much less than M , the probability 
of undercounting due to imperfect splitting will be small. 
However, when N is comparable to M , we must account 
for splitting errors. Using a standard urn model |18LI21j . 
we can calculate the conditional probability of detecting 
n photons given that the incoming field has N photons 
in it. This probability is given by the expression 



Pdct{n\N) = 
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where rj is the single-photon detection efficiency. It is not 
difficult to verify that in the limit M » 1 and M > 
detection errors due to imperfect splitting are negligible 
and that Pdctin\N) -> (^) (1 - r;)^""?]", which is the ex- 
pected distribution for an ideal photon-number-resolving 
detector subject to a loss of 1 — 77 per photon. We note 
that detector dark counts are a source of overcounting 
errors. However, in the present analysis, we do not in- 
clude the effect of detector dark counts, since often they 
can be rendered neglible through detector gating or by 
use of detectors with intrinsically low dark count rates 
(commercially available APD dark count rates can be as 
low as 5 dark counts per second). 

The response of the detection system is shown in Fig. [2] 
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FIG. 2: Regions where there is a significant probability of 
detecting n photons conditioned on A'^ photons being incident 
on the TMD for three different single-photon detection effi- 
ciencies (7; = 1.0, 0.66, 0.33). The mean number of photons 
detected given N photons input into the detection system are 
plotted using the solid {rj — 1), dotted {rj = 0.66) and dashed 
{ri = 0.33) lines. The shaded regions are within one standard 
deviation of the means. 



for a system comprised of five beam splitters (M=32) for 
three different scenarios. For no loss in the TMD (77 = I) 
and few enough incident photons, the detector will accu- 
rately measure the number of photons in the field, i.e. the 
initial slope of the mean of the probability distribution 
is linear with unity slope. Once the number of incident 
photons increases and becomes comparable to the max- 
imum number of photons M the system can count, un- 
dercounting will become very likely and the detection 
system will saturate. As the loss in the system increases 
(77 decreases), two effects are seen: (I) the initial slope 
decreases and is given roughly by 77 and (2) the error 
in the estimate increases as seen by the increase in the 
vertical spread of the distribution in the plot. 

As we already mentioned, in a heralding system any 
error in the detection of the number of photons in the 
heralding field will directly lead to an uncertainty in the 
number of the photons in the heralded field. If this un- 
certainty, as determined by the variance of the number 
photons estimated to be in the heralding field, becomes 
larger than the mean number of photons in the heralded 
field, the state can, in principle, be produced by classical 
means. However, since the source of the photon pairs 
is an OPA, we can use this a priori information about 
its photon emission statistics to improve our estimation 
through the use of Bayesian theory. 

The probability that an OPA characterized by gain 
g will emit k pairs of photons is given by PoPA{k) — 
\{k, k\ip)\^ , where \k,k) is the state containing k pho- 
tons in both the signal and idler fields and is the 
state emitted by the OPA, which is given by exp{ga\al — 
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FIG. 3: Comparison of the photon-number statistics for a 
typical heralded field with n™' = 5 (in blue) and a Poisson- 
distributed field with a mean photon number of 5 (in gray). 
The heralded field, for the case in which r) — 0.66 and 5=1, 
is conditioned on the detection of 4 photons in the idler field 
by the TMD. 



(70^0,3)10,0) [22]. As a notational device, we will identify 
the signal field with the heralded field and the idler field 
with the heralding field. It can be shown that the result- 
ing photon-number distribution for either the signal or 
idler field obey Bose- Einstein statistics and is given by 



PoPA(fc) — (1 ^ tanh'^g) tanh'^'^g 



(2) 



We can calculate the resulting mixed state of the signal 
field given that rii photons are detected in the idler field. 
The density matrix is given by 



Psig("i) = X! ^sig('^sl'^i) |?^s)(«s|, 



(3) 



where jrig) is a Fock state of the signal field with rig pho- 
tons in it. The function Psigins\ni) is the probability that 
the signal field has Ug photons in it given that ni pho- 
tons are detected in the idler field. According to Bayes' 
theorem, it is given by the expression 



-Psig(»^s|"i) 



■Pdct("-i|"'s)-PoPA('^s) 
Lfclo Pdet{ni\k) PoPA(fc) ' 



(4) 



We estimate the number of photons in the signal field 
using maximum-likelihood estimation. The maximum- 
likelihood estimate of the number of photons in the signal 
field given rij photons detected in the idler field is found 
using 



n'^\ni) = argmaxPsig(ns|ni), 



(5) 



where the argmax function returns the value of rtg G 
{rii, Hi + I, ... , 00} that maximizes Psig('^s|"-i)- 

In Fig. |3] we show the photon-number statistics for 
a heralded field estimated to have n™' = 5 photons 
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FIG. 4: Uncertainty in the number of photons in a her- 
alded field when ni photons are detected in the heralding 
field normalized to the uncertainty in a comparable Poisson- 
distributed field with mean number of photons n™' for vari- 
ous detection efficiencies {rj — 0.33,0.66, 1.0) and OPA gains 
(5 = 0.75,1). 



in it and compare it to a classical coherent state that 
has Poisson statistics with a mean number of photons 
equal to 5. From the figure it is clear that the photon- 
number distribution of the heralded state is much nar- 
rower than the distribution for the comparable classical 
coherent state. The width of the heralded field distri- 
bution is exactly the error in the maximum-likelihood 
estimate and is found using the mean squared error, 
i.e. (A<i)'. = EnMs-nffPsisin^lm). FigureHcom- 
pares the ratio of the spreads of the heralded fields to 
their classical counterparts for different numbers of pho- 
tons in the heralding field for various values of the single- 
photon detection efficiency rj and the OPA gain g. When 
this ratio is equal to unity, the heralded field has as much 
uncertainty in its number of photons as a Poissonian field. 
The smaller the ratio is, the closer the heralded field is to 
a Fock state. The quality of the heralded state degrades 
as the loss increases and the gain of the OPA increases, 
as can be seen by noting that the ratio increases in these 
cases. 

To gain further insight into the behavior of the her- 
alded state, we consider the mean number of photons in 
the heralded field. Although the mean number of pho- 
tons cannot be an estimator of the number of photons in 
the heralded field, since it can take non-integer values, it 
is useful parameter that characterizes the field and can 
be analytically calculated. The mean number of photons 
in the signal field conditioned on ni photons being de- 
tected in the idler field can be expressed in the following 
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remarkably simple form (see Appendix) 

Similarly, the conditional variance of the number of pho- 
tons in the signal field can be expressed by 

(An,)^^^^(4-a,), (7) 

where in both equations we have introduced the quan- 
tity = {l — tanh^(7 [(1 — rj) + ri{n — j) /M] } for no- 
tational convenience. An ideal system would produce 
Fock states with definite photon number so that (Arig)^. 
would vanish. We note that aj cannot vanish since that 
requires rj to vanish, which is not a physically meaningful 
situation. Thus, (Ang)^. goes to zero only in the limit 
in which aj I. Also in this limit the conditional mean 
{ns)ni ^ n-i, which is consistent with the expected be- 
havior of an ideal heralding system. 

Two cases exist in which aj approaches unity. The 
first case occurs in the low gain limit. Irrespective of the 
system loss t], aj can be made arbitrarily close to one for 
small enough gain since aj < 1/(1 — tanh^t;) = l + 0{g^). 
Intuitively, this limit is easily understood. When the 
gain is small enough, it is much more likely for the 
source to produce rig photon pairs and the detector to 
count rii of them than it is for the source to produce 
Us + I pairs and the detector to still detect only rtj. 
More quantitatively, there exists a gain such that the ra- 
tio r = Pdot(n-i|".s)-PoPA('^s)/(^'dct(".i|'^s + 1)Popa("-s + 
1)) 3> 1. This statement is easily verified by noting that 
r c>c l/tanh^(7 and can be made large for small enough 
g. The second case in which aj goes to one occurs when 
M ^ n and i] — 1. In this limit the TMD behaves like an 
ideal photon number resolving detector and the system 
achieves perfect heralding. 

Although these limiting cases are useful in providing 
insight into the behavior of the system, they are often 
not viable experimentally. If the gain is too low, the rate 
at which the heralded states are produced becomes im- 
practically small. Furthermore in many cases M cannot 
be too large, if for no other reason than that the loss in- 
curred by adding splitters can become excessive. To treat 
the problem more generally, we calculate Mandel's Q pa- 
rameter defined hy Q = ((Ans)^. - (ns)„i) /(ris)™; [I]- 
Mandel's Q parameter is useful for characterizing the na- 
ture of the photon-number distribution of a field. The 
Q parameter is bounded below by —1. In the limit in 
which Q = — 1, the heralded field is in a Fock state. For 
Q = 0, the field has Poissonian photon-number statistics; 
for positive and negative Q the field has super- and sub- 
Poissonian statistics respectively. We use this parameter 
as a merit function describing the quality of the heralded 
states produced using the technique. 
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FIG. 5: Mandel's Q parameter versus the gain of the 
OPA and the single-photon detection efficiency of the time- 
multiplexed detector (TMD) for the case where n™' = 5, i.e., 
the heralded field is a single-mode field estimated to contain 
5 photons. The colored regions are labeled with the num- 
ber of detected idler photons needed to herald the presence 
of 5 photons in the signal field. The white contours are lines 
of constant P{g,ri;n^^), the probability of creating the her- 
alded state. The TMD is assumed to contain 5 beam spUtters 
(M=32) 

In practice, it is useful to know the TMD and OPA 
parameters necessary to produce a state that contains 
an estimated n™' photons and a certain value of the Q 
parameter. To this end, we calculate Q at constant ' 
versus the parameters rj and g. Figure [5] shows the result 
of the calculation for a heralding detector comprised of 
five beam splitters for the case when ' = 5. To calcu- 
late Q for constant n™' it is necessary to invert Eq ^ 
to find the appropriate value of rij to be substituted into 
Eqs ([6| and (|7| . For certain values of the gain and detec- 
tor efficiencies, Eq ^ cannot be inverted, i.e. there is no 
integer number of detected photons ni that correspond 
to a signal field in which n™' = 5. In these instances, 
a signal field estimated to contain five photons is not 
possible to produce and is marked by the gray regions 
in the plot. Also, the gray bands divide the space into 
regions corresponding to the number of detected idler 
photons that herald the presence of 5 photons in the sig- 
nal field. The figure confirms that the losses as well as 
the OPA gain should be kept as small as possible. In 
addition, the efficiency with which the heralded states 
are created is important. The white contour lines over- 
laid on the plot show the probability with which the her- 
alded state is successfully created, which is denoted by 
P{g, rj; n™'). This probability is calculated from the sum 
SfcLo -Pdet('^ilfc) ^'oPA(fc) where is obtained from the 
inversion of Eq. (|5| for different values of g and rj. The 
discrete jumps seen in the contour lines are due to the 
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FIG. 6: Mandel's Q parameter versus gain of the OPA and 
the single-photon detection efficiency of the TMD for the case 
where n™' — 5 and = 5, i.e. the signal field is estimated 
to contain 5 photons spread between 5 modes. The colored 
regions are labeled with the number of detected idler pho- 
tons needed to herald the presence of 5 photons in the signal 
field. The white contours are lines of constant P{g,r];n^^), 
the probability of creating the heralding state. 

fact the Hi can take only integer values. We see that it 
is possible to create nearly ideal Fock states even with 
significant detector loss; however, the efficiency of the 
creation of the heralded states is sacrificed. 

Up to this point, we have considered only single-mode 
states produced by an OPA. Although the single-mode 
limit can be reached, e.g. by both spatially and fre- 
quency filtering [53] , the photons emitted by a free-space 
OPA are generally multimodal in nature and can be ap- 
proximately described by a multimode thermal distribu- 
tion [23] . The probability that the OPA will emit k pairs 
of photons distributed among any modes is given by 

PoPAAk) = ^ ^ 1 ^) (1 - tanh^g)'^ tanh^'^s, (8) 

where g is the parameter characterizing the gain per 
mode of the OPA. We have assumed that each mode 
is equally likely to be occupied. The multimode ther- 
mal distribution has a variance of (fc) (1 -I- Since 
when (k) is held constant, the variance decreases with 
increasing number of modes, it is reasonable to expect 
that we may be able to create multimode fields contain- 
ing near-definite numbers of photons under more relaxed 
conditions than by restricting ourselves to single-mode 
states. 

In Fig. |6]we plot the Q parameter under the same con- 
ditions as in Fig. [5] but now assuming that we are inter- 
ested in the multimode output of the OPA (see Appendix 
for calculation details). We take = 5 as an example. 



which is typical value for a free space OPA [2S|. When 
comparing Fig. [5] and Fig. (e] we see that the multimode 
case is more tolerant of loss in the heralding detector 
than the single mode case. For example, if we require a 
heralded state preparation rate of 5%, the heralded out- 
put of the OPA is sub-Poissonian so long as rj > 0.11, 
but in the single-mode case we have the more stringent 
requirement that rj > 0.36. 

In conclusion, we have shown that it is possible to cre- 
ate multiphoton heralded states with highly non-classical 
photon-number distributions using a TMD as the herald- 
ing detector. We calculated Mandel's Q parameter as a 
function of the detector loss and the OPA gain to charac- 
terize the quality of the heralded state. Even in the pres- 
ence of significant loss, we found that it is possible to pro- 
duce states with highly sub-Poissonian photon-number 
distributions at reasonable efficiencies. In addition, for 
situations where a single-mode output is not required, we 
showed that the full multimode output of an OPA can 
provide less uncertainty in the number of photons in the 
field than the single-mode output. Since the TMDs in- 
volve only commercially available APDs and fibers, the 
technique we have presented here can be a particularly 
simple way to create states of light containing nearly def- 
inite numbers of photons. 

The authors would like to thank Christine Silberhorn 
and Gerd Leuchs for helpful discussions. This work 
was supported by the US Army Research Office through 
a MURI grant. KWC thanks the support from the 
Croucher Foundation. 

APPENDIX 

The derivation of the conditional mean and variance of 
the number of photons in the signal field conditioned on 
rij photons being detected in the idler field is described 
in this appendix. 

To calculate the conditional mean and variance, we 
need to evaluate the sums defined by 

oo 

Si=J2 niPdet{ni\ns)PoPAAns), ^ = 0, 1, 2, (A.9) 

n^=Q 

where Pdct('T-il'T-s) and PoPA,/i('T-s) are the probability dis- 
tributions given in the text in Eq. ([T]) and Eq. (|8| respec- 
tively. 

The simplest case is that of / = 0. Performing the 
summation over Us by recognizing the negative binomial 
series, we find 

so=Af:h)jtiy, (A.io) 

where we have defined A = (*^) (^'//(rysinh^.g))'' and 
X = M / (rj siiih^ g) + M — ni for simplicity. Next, 
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using the identity for Gauss hypergeometric series, 
E;=o (7)(-1)V(^ + j) = + where B(a,6) = 

T{a)T{b)/T{a + b) is the beta function, we express as 



where dj^f^ is shorthand for the (/i — l)th-order partial 
derivative with respect to x. 

We follow an analogous procedure for the case when 
I — 1,2, where the primary difference is that the summa- 
tions over Ug are now derivatives of the negative binomial 
series. The resulting sum 5*1 is given by 

= [{x + b) B{x, 1 + m)] , (A. 12) 

with b = Af (1 — 7^)/?/ + rij. The sum S2 is similarly given 

by 



^'-^ (M+1)! ^ 



{x + bf B{x,l + ni) 



(x + b)c 



B{x, 1 + rii) 



(A.13) 



where c — M / {rj tanh^ g) is defined for simplicity. 

The conditional mean and variance can be expressed 
in terms of these sums as (ns)„i — Si/Sq and (Ang)^. = 
82/80 — (S'l/S'o)^. Substituting into these equations and 
simplifying, we find the general form for the conditional 
mean and variance to be 

(ns)„i = -^-cgf_,{x) (A. 14) 

(Aris)^, = {c + c^d,)gf,{x), (A.15) 

where 5^(x) = di^'^ B{x, 1 + n{) / di^^^'^ B{x, l + m). By 
setting /i = 1 it is not difficult to recover the single- 
mode results of Eq. ^ and Eq. (|7| shown in the text. 
Furthermore, the single-mode and multimode expressions 
for Q can be easily found from these results. 
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